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Abstract 

In this paper, a new method to determine the complete set of dynam- 
ically balanced planar four-bar mechanims is presented. Using complex 
variables to model the kinematics of the mechanism, the dynamic balanc- 
ing constraints are written as algebraic equations over complex variables 
and joint angular velocities. After elimination of the joint angular ve- 
locity variables, the problem is formulated as a problem of factorization 
of Laurent polynomials. Using toric polynomial division, necessary and 
sufficient conditions for dynamic balancing of planar four-bar mechanisms 
are derived. 



1 Introduction 

Statically and dynamically balanced mechanisms are highly desirable for many 
engineering applications since they do not apply any forces or moments at their 
base, for arbitrary motion trajectories. This concept is used in mechanism 
design in order to reduce fatigue, vibrations and wear. Static and dynamic bal- 
ancing can also be used in more advanced applications such as, for instance, the 
design of more efficient flight simulators [5], or in the design of compensation 
mechanisms for telescopes. Additionally, dynamic balancing is very attractive 
for space applications since the reaction forces induced at the base of space ma- 
nipulators or mechanisms are one of the reasons why the latter are constrained 
to move very slowly [16]. 
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under the SFB grant F1303. B. Moore was also supported by the Fonds Quebecois de la 
Recherche sur la Nature et les Technologies (FQRNT). C. Gosselin was partially supported 
by the Natural Sciences and Engineering Research Council of Canada (NSERC). 
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Several approaches can be used to balance mechanisms (see for instance 
[2, 1]). In general, complete balancing requires the integration of additional me- 
chanical components in the design of a mechanism, such as counterweights and 
counterrotations [4]. However, for some simple architectures, it is sometimes 
possible to design dynamically balanced mechanisms by an appropriate choice 
of the design parameters without introducing additional linkages or counterrota- 
tions [9]. Even though these mechanisms do not include counterrotations, they 
satisfy all conditions for balancing, namely, their centre of mass remains sta- 
tionary (static balancing) and their total angular momentum vanishes (dynamic 
balancing), for arbitrary trajectories. 

Although families of dynamically balanced four-bar mechanisms were pre- 
sented in [9], no proof was given on the existence of other possible solutions. 
In this paper, the aim is to derive all possible sets of design parameters for 
which a planar four-bar mechanism is dynamically balanced. The problem was 
first addressed by Berkof and Lowen [3] , who provided conditions for static bal- 
ancing in terms of the design parameters when the geometric parameters are 
sufficiently generic. A non-generic solution was then found in [8]. In [11], a 
complete list of statically balanced planar four-bar mechanisms was given. The 
problem of dynamic balancing was also addressed in [14], where special cases 
that do not require external counterrotations were first revealed. In the latter 
reference, it can be observed that the dynamic balancing problem leads to a 
rather complicated system of algebraic equations. 

A generic approach to investigate such parametric polynomial systems has 
been recently proposed by Lazard and Rouillier [10]. In [11], it was shown that 
the problem can be simplified if one uses complex variables to model the angles 
in the configuration space, together with some results in algebraic geometry 
(Ostrowski's theorem [12]). However, the application of the above techniques to 
the problem of dynamic balancing leads to a cumbersome case by case analysis. 
In order to simplify this analysis, a concept of division for Laurent polynomials 
is introduced here. The method is related to [15], who use similar ideas for 
factoring polynomials by taking advantage of the special shape of their Newton 
polyhedra. 

This paper is organized as follows. First, we describe how to derive the 
kinematic, static and dynamic equations using complex number representations. 
This leads to a set of algebraic equations in two complex variables where the 
coefficients are expressions in terms of the design parameters. This section 
can be skipped if one is interested only in the mathematical aspects. Then, 
we introduce the concept of division for Laurent polynomials and provide an 
algorithm for computing this division. Then, we use Laurent polynomial division 
to eliminate the dependent variables describing the configuration space. Finally, 
we solve the system in the remaining design parameters and give an overview 
of the various subcases. 

2 Model 

2.1 Representation of planar four-bar mechanisms 

A planar four-bar mechanism is shown in Figure 1. It consists of four links: 
the base of length d which is fixed, and three moveable links of length h,h,h 
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Figure 1: Four-bar mechanism. 



Type 




Parameters 


Kinematic 


Length 


h,h, h, d 


Static 


Mass 
Centre of mass 


mi,TO2,m 3 
n,ipi,r2,i>2,r3,ip3 


Dynamic 


Inertia 





Table 1: Design parameters for the planar four-bar mechanisms. 



respectively. We assume that all link lengths are strictly positive. Since the 
base is fixed, the mass properties of the base has no influence on the equations 
and will therefore be ignored. Each of the three moveable links has a mass m i} 
a centre of mass whose position is defined by and ipi and a moment of inertia 
Ii. The design of planar four-bar mechanisms consists in choosing the 16 design 
parameters (Table 1). 

The links are connected by revolute joints rotating about axes pointing in a 
direction orthogonal to the plane of motion. The joint angles are specified using 
the time variables 9i(t),02(t) and 9 3 (t) as shown in Figure (1). Since the mech- 
anism has only one degree of freedom, there is a relationship between these joint 
angles, which will be described below. The kinematics of planar mechanisms 
can be conveniently represented in the complex plane, using complex numbers 
to describe the mechanism's configuration (Figure 2) and the location of the 
centre of mass (Figure 3). Referring to Figures 2 and 3, let zi,Z2,Z3 be time 
dependent unit complex numbers and pi, P2, P3 unit complex numbers depend- 
ing on the design parameters (actually only on tpi, tp2, ^3). The orientation of 
Pi is specified relative to Zj, i.e., it is attached to z^ and moves with it. If p, 
coincides with z i; then p» = 1. 
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2.2 Kinematic model 

The dependency between the different joint angles is described by the following 
closure constraint: 

z 3 = Gizi + G 2 z 2 + G 3 (1) 

where G 1 ,G 2 ,G 3 eR with Gi = G 2 = ff , G 3 = Taking the time deriva- 
tive of equation (1), we get a relationship between the joint angular velocities 
6*i, #2 and 63, namely: 

z 3 03 - Gizifli + G 2 z 2 2 (2) 

Since Z3 is a unit complex number, Z3Z3 = z 3 z 3 1 = 1 and therefore we 
obtain the following geometric constraint: 

G = (Gizi + G 2 z 2 + G 3 ) (Gxz- 1 + G 2 z 2 " 1 + G 3 ) - 1 = 0. (3) 

The time derivative of the geometric constraint (3) can be written as a linear 
combination of the joint angular velocities: 

i{KA + K 2 6 2 ) = (4) 

where 

K x = GiG 2 (ziZ2 1 - zr 1 z 2 ) + GiG 3 (zi - O (5) 

if 2 = G 1 G 2 (z- 1 z 2 - z^- 1 ) + G 2 G 3 (z 2 - z- 1 ) (6) 

It is noted that since K\ and K 2 are purely imaginary, only one constraint 
equation is obtained, over the real set. 
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Figure 4: Unit vectors representation 



2.3 Position of the centre of mass 

Let M be the total mass of the mechanism (M = mi +1112 + 1713). The centre 
of mass C is: 



C = (Fizi + F 2 z 2 + F 3 ) 



where Fi,F 2 ,F 3 e C: 

Fi 
F 2 
F 3 



miripi + m^li + Gim 3 r 3 p 3 
TO2?"2P2 + G 2 m 3 r 3 p 3 
m 2 d + G 3 m 3 r 3 p 3 . 



These equations were derived in [ll] 1 . 



(7) 



(8) 
(9) 
(10) 



2.4 Angular momentum of the mechanism 

Since the mechanism is planar, the contribution of body i to the angular mo- 
mentum is a scalar and can be given in the following form: 



H, 



r i/o,~ir i/o) +1$. 



(11) 



where T^/q and ~~f' \/ are respectively the position and the velocity of the centre 
of mass of body i with respect to a given inertial frame, Ii denotes the moment of 
inertia of body i with respect to its centre of mass and (*, *) is the scalar product 
of planar vectors. The total angular momentum H of the system is given by 
the sum of the angular momentum of the links, (i.e. H — H 1 + H 2 + H 3 ). 



1 In this paper, indices 2 and 3 for the links have been permuted. 
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The angular momentum of the first body with respect to the inertial frame 

is: 



Hi = (ripizi,-imi(iripizi0i)) + h9i = (ripizi, minpizifl'i) + Ji 6*i = Ji6\ 

(12) 

The contribution of the second body to the angular momentum is given by: 



H 2 = {(d + r 2 p 2 Z 2 ) ,m 2 r 2 p 2 Z 2 2 )+I 2 6 2 



/ p 2 z 2 + p 2 x z 2 1 
m 2 dr 2 I 



For the third body, we get 

H 3 = ( (Zizi + r 3 p 3 z 3 ) , m 3 ( ZiZi#i + r 3 p 3 z 3 6> 3 ) ) + I 3 6 3 



2 
(13) 



(14) 



Substituting equations (1) and (2) into equation (14), we can eliminate z 3 # 3 
and 9 3 and obtain en expression in terms of zi, z 2 , 9\,9 2 only. The total angular 
momentum H of the mechanism is then given by: 

H = Hi + H 2 + H 3 = K 3 9\ + Kj 2 (15) 

where K 3 and K4 are written as: 

K 3 = aizi + a^ 1 + foiziz^ 1 + & 2 z 1 ~ 1 z 2 + c 

if 4 = UlZ 2 + U2Z2 -1 + WiZiZt" 1 + W 2 Z^ 1 Z 2 + W 

where constants a\,a 2 ,b\,b 2 , c, u\, u 2 , v\,v 2 and u; can be obtained from equa- 
tions (12, 13, 14, 15). 

2.5 Static and dynamic balancing 

In our settings, a mechanism is said to be statically balanced if the centre of 
mass of the mechanism remains stationary for infinitely many configurations (i.e. 
infinitely many choices of the joint angles). From equation (7), this condition 
can be formulated as: 

F = Fizi +F 2 z 2 - C = (16) 

where C — CM — F 3 is a constant. A mechanism is said to be dynamically 
balanced if the total angular momentum remains constant for any motion of 
the mechanism. In other words, the reaction forces and torques at its base 
induced by its motion are identically equal to 0, at all times. Clearly, the 
mechanism must be statically balanced and the angular momentum must be 
constant. Since the centre of mass is fixed and we want static balancing for all 
possible configurations and all joint angular velocities, the angular momentum 
should therefore be 0, i.e., 



H = K 



3 o\ 



4^2 







(17) 



Therefore, equations (3, 4, 16, 17) have to be satisfied. Among these four 
equations, only two (equation 4 and 17) depend (linearly) on the joint angular 
velocities and they can be rewritten in the following form: 



(18) 
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If the rank of the matrix A — 



is 2 and since the system is homo- 



K x K 2 

K 3 Ki 

geneous, then the only solution is Q\ = 9 2 = 0. In other words, the mechanism 
is not moving. Therefore we must have: 



K := det(A) = K X K A - K 2 K 3 = 



(19) 



Using this equation together with the geometric and static balancing con- 
straints, the joint angular velocities (6i, 2 ) are eliminated and a system of three 
algebraic equations (equation 3,16,19) in terms of the unit complex variable 
Zi,z 2 is obtained. In section 4, all possible design parameters satisfying this 
system of equations for infinitely many configurations will be derived. However, 
we first need to introduce concepts and tools from toric geometry in order to 
solve these equations. This is the subject of section 3. 



3 Toric geometry 

Definition 1 A Laurent polynomial g over a ring R is a formal sum of mono- 
mials x a := Cq.x" 1 ^^ 2 , where x = (xi,x 2 ) is a fixed pair of variables, and 
a G 1? , c a G R. Its support is the set of all a £ Z 2 with non-zero coefficients 
c a . Its Newton polygon is the convex hull of the support in M 2 . The Laurent 
polynomials form a ring, namely R[xi,x 2l (xix 2 ) -1 ]. 

Definition 2 The Minkowski sum of two convex sets A and B cR™ is defined 
as: 

A + B = {a + b\aeAAbeB} (20) 
Note that A + B is also a convex set. 

Theorem 1 Assume that R does not have zero divisors 
polynomials, then: 

U(FG) = 11(F) + 11(G) 

We refer to Ostrowski, [12, 13] for a proof. 

Remark 1 The assumption that R has no zero divisor can be replaced by the 
weaker assumption that the corner coefficients of G, i.e., the coefficients at the 
vertices ofH(G), are not zero divisors. 

In order to find out whether a given polynomial G divides another given 
polynomial F, we introduce Laurent polynomial division. 

Definition 3 Assume that G is a Laurent polynomial such that its corner co- 
efficients are no zero divisors. A finite subset T of 1? is called a remainder 
support set with respect to G iff no multiple of G, except zero, has support 
contained in V. 

Definition 4 Let G be a Laurent polynomial such that its corner coefficients 
are invertible in R. Let F be an arbitrary Laurent polynomial. Then (Q, R) is 
a quotient remainder pair for (F, G) iff the following conditions are fulfilled. 

a) F = QG + R. 

b) The support of R is contained in 11(F). 

c) The support of R is a remainder support set with respect to G. 



. IfF,G are two Laurent 
(21) 
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Algorithm 1 Toric Polynomial Division Algorithm 



1: Input: F, G, such that G has invertible corner coefficients. 
2: while F^Odo 

3: Select a linear functional h : M 2 — > K, (a;, y) (as + fly) such that a//3 
is irrational. 

Compute the point / e Support(F) that maximizes h. 
Compute the point g e Support (G) that maximizes h. 
if 11(G) +f-gc U(F) then 

Q:=Q + M;F:=F- MG; 
else 

M := Coef(F,/)a^V 2 ; 
i? i? + M; F := F - M; 
end if 
end while 
Output: Q, i?. 



Quotient remainder pairs are not unique. Here is a nondeterministic algo- 
rithm that computes quotient remainder pairs. 

Theorem 2 Algorithm 1 is correct. 

Proof 1 The Newton polygon of F becomes smaller in each while loop, hence 
it is clear that Algorithm 1 terminates. Also, any monomial which is added to 
R is contained in 13(F), hence it follows that R fulfills (b) in Definition 4- 

No step in the algorithm changes the value of F + QG + R. Initially, this 
value is the given polynomial F, and in the end, this value is equal to QG + R. 
This shows that (a) in Definition 4 is fulfilled. 

In order to prove (c) in Definition 4, we claim that the following is true 
throughout the execution of the algorithm: if H is any Laurent polynomial such 
that GH has support in n(F) U Support(F), then the coefficients of GH at the 
exponent vectors in Support(F) are zero. 

Initially, Support (R) is empty and the claim is trivially true. If the claim 
is true before step 8, then it is also true after step 8, because this step does not 
change R and does not increase the Newton polygon of F. 

Assume that for a certain Laurent polynomial H , the claim is true before 
step 11 and false after step 11. Then it follows that the coefficient of GH at 
f is not zero, because this is the only exponent vector which is new in R. The 
support of GH is also contained in the Newton polygon of F before step 11, 
hence f is the unique vector in Support (GH) where h reaches maximal value. 
Because g is the unique vector in Support (G) where h reaches a maximal value, 
it follows that (/ - g) G Support (H). Then 11(G) + f - g C T1(GH) as a 
consequence of Theorem 1. But this implies that the if condition in step 6 is 
fulfilled for G and F before step 11, and therefore step 11 is not reached for such 
values of F and G. 

It follows that the claim is true throughout the execution of Algorithm 1. In 
particular, it is true at the end, which shows that (c) in Definition 4 holds. 

Example 1 Let F = c 02 y 2 + ciixy + c iy + c w x + c 00 and G = d iy + d 10 x + d 00 . 
The result of the polynomial division algorithm is shown in table (2). Therefore, 
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Step 


f 


g 


Computation 


7: 

8: 


(0,2) 


(0,1) 


dm y 
^ dm a 

F = (en - a^ui) xy + (coi - ^20) y + c 10 a; + c 00 


10: 
11: 


(1,1) 


(0,1) 


M = ( c " - ^) x v 

R — ( r c 02 d 10 \ r7 . 

n — \ Ci1 d m ) x y 

F = (coi - y + c 10 x + c 00 


7: 

8: 


(0,1) 


(0,1) 


Q = i^y + d^ 


10: 
11: 


(1,0) 


(0,1) 




10: 
11: 


(0,0) 


(0,1) 


R = ( c " - «W + ( c io - 4ft 1 ) * + ( c «o - ) 
F = 



Table 2: Example: Toric polynomial division algorithm 



F is divisible by G if and only if R — 0, or in other words if all coefficients of 
R are zero: 

£02^10 Adw Ad$Q 

en = C10 - — — = coo - — — = (22) 

"01 «01 "01 



4 Balancing 

4.1 Problem description 

The problem addressed in this paper can be stated as follows: find all possible 
design parameters such that there exists a valid non-constant trajectory of the 
planar four-bar mechanism for which the mechanism is dynamically balanced 
(i.e.: it is statically balanced and the angular momentum of the system is 0). 
Formally, let K = {(zi,z 2 € C 2 | G(zi,z 2 ) = 0)} be an infinite set representing 
a valid non-constant trajectory. For this trajectory, we want the mechanism to 
be statically balanced, i.e.: 

V (zi;Z2)eK G(zi,z 2 ) = => F(zi,z 2 ) = (23) 

and the angular momentum to be 0: 

V (zi;Z2)eK G(zi,z 2 ) = => K(z 1 ,z 2 ) = (24) 

where G, F and K are defined in equations (3, 16, 19). Their Newton polygons 
are shown in figure 5. Using the following theorem, we can reformulate this 
problem as a factorization problem of Laurent polynomials. 
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P.-2) 



n(G) 11(F) U(K) 

Figure 5: Newton polygons of G, F and K. 

Theorem 3 Let G be an irreducible Laurent polynomial. Let K be a Laurent 
polynomial (not necessarily irreducible). Let S C C* 2 such that G has infinitely 
many zeros in S . The following are equivalent: 

1. V( Zl ,z 2 ) S S,G(zi,z 2 ) = 0=S> K(zi,z 2 ) =0 

2. 3 Laurent polynomial L such that K = G ■ L 

A proof of this theorem can be found in [11]. 

4.2 Static balancing 

Assume G(zi,Z2) is irreducible and using theorem (3) we are looking for a 
Laurent polynomials L(zi,z 2 ) such that 

F(z 1 ,z 2 ) = G(z 1 ,z 2 )L(z 1 ,z 2 ) (25) 

If the geometric constraint G is not irreducible, we consider all possible de- 
composition of G into irreducible components (see table 3) . Every such decom- 
position imposes constraints on the kinematic parameters /i,Z 2 ^3 and d (table 
4). For a given decomposition, one component corresponds to a kinematic mode 
of the mechanical system. For each of these decompositions and components, 
we can apply theorem (3). Using this approach, necessary and sufficient condi- 
tions for the static balancing of planar four-bar mechanisms can be obtained as 
shown in [11]. These conditions are described in table 6. 

4.3 Dynamic balancing 

The same approach as in the static balancing case could be used to find all suffi- 
cient and necessary condition for the dynamic balancing (equation 24). However, 
due to the complexity of the Newton polygon of K, the approach would lead to 
a cumbersome case by case analysis, making it unpractical and prone to error. 
Using the toric polynomial division algorithm 1, the same result can be obtained 
in a semi-automatic way (using symbolic computation tools) without this case 
by case analysis. In our case, the computations were performed using a Maple 
implementation of the Toric Polynomial Division algorithm and the solutions 
are presented below. 
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^ ^ ■ □ 

(I) (II) 



(in) (iv) 




(V) 

Table 3: Possible decomposition of G based on Newton polytopes and Minkowski 
sums. 



Case 


Kinematic constraint 


Mode A 


Mode B 


Mode C 


II 


h — h,h — d, h 7^ ^3 


Zl 7^ z 2 


/ 7 

Zl = z 2 




III 


W 

^i — d, I2 — I3, l\ 7^ h 


zi 7^ 1 


/ 

zi = 1 




IV 


h = h,h = d, h 7^ I2 


z 2 7^ 1 


J 

7,2 = 1 




V 


h=h = h = d 


7.2 = —i 


_y 

zi = 1 


Zl = z 2 



Table 4: Kinematic mode. 
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4.3.1 Irreducible case 

Figure 6 illustrates the main steps of the polynomial division algorithm when the 
geometric constraint G is irreducible. The algorithm gives a set of constraints in 
terms of the design parameters which can be combined with the static balancing 
constraints. Among these constraints, we obtain 

w = lp 3 + lp 2 = (26) 

with J2 = m-2,r\ + I2 and J3 = m?,r\ + I3. Therefore I2 = h = r 2 = ^3 = 
which is physically not possible. Therefore, if G is irreducible, a planar four-bar 
mechanism cannot be dynamically balanced. 

4.3.2 Reducible case II, Mode A 

Here we have I3 = d and l\ = h =■ I- This is one case where we get solutions 
that are physically realizable. In order to decribc the solutions, wc introduce 
another set of parameters, namely 

<7i := npiini, Ji := I { + mirf,i = 1, 2, 3. 

The variables Ii , r j , Pi can then be eliminated easily. The balancing problems 
have an additive structure: if a balanced mechanism picks up weights at each of 
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h = 1 


h = l 


Z 3 = 4 


d = 4 


mi = mi 


m 2 = ^ 


™3 = ^f 1 




n = h 




^3 = 1 




Pi = -1 


= -1 


P3 = 1 




j 3m! 









Table 5: Example of a dynamically balanced mechanism in case 2 A 



its three bars in a balanced way, then the composed mechanism is also dynam- 
ically balanced. The variables are now chosen so that the balancing conditions 
become linear in raj, qt, Ji, i = 1, 2, 3. Here they are. 

Qi = ^<73 - lm 3 , q 2 = -^3, (27) 

= — ^ — 93 - Ji - I m 3 , J 2 = — - — q 3 -J3- (28) 

A first consequence is that (71,(72, 93 must be real. The parameters fulfill also 
the inequality constraints 

m 4 > 0, J 4 m 4 - | ft | 2 > (29) 

for i = 1,2,3. In particular, J x and J 2 must be positive. By equation (28), we 
get an upper bound for m 3 , which must be larger than the lower bound from 
(29) (i = 3). This yields 

(dq 3 - J 3 )(dJ 3 -l 2 q 3 ) >0. (30) 

It follows that J 3 is contained in the open interval {dq 3l l ^q 3 ). (Note that 

d ~^ q 3 > as a consequence of (28), that is why we know which of the two 
interval boundaries is bigger.) Then q 3 > and d > I follows. From J 2 > and 
(28), we get 

12 T d2 - 12 

1 93 < J3 < (31) 

from which d > \/~2l follows. 

Conversely, if d > V%1, then we can choose q 3 > arbitrarily and J3 subject 
to (31), and m 3 between the upper and lower bound for m 3 derived above. Then 
(28) determines Ji and J2, which will then be positive. Then (27) determines q\ 
and (72, and finally mi and m.2 can be chosen so that inequality (29) is fulfilled. 
A possible solution is given in table (5). 

4.3.3 Reducible case II, Mode B 

In this kinematic mode, the mechanism is a parallelogram with z\ — Re- 
placing zi by z 2 , the constraint K = becomes: 

~~~T2~ (z2-l)(z2+l) [(01 + ?ii)z 2 + (61 + b 2 + c + vi + v 2 + w)z 2 +u 2 + a 2 ] = 

2 

" (32) 

Therefore, the angular momentum vanishes if and only if one of these factors 
vanishes. The first two factors (i.e.: z 2 — 1 and z 2 + 1) correspond to uncertainty 
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configurations in which it is possible to switch between mode A and mode B 
(i.e.: where we can pass from one mode to the other) and are valid for only one 
configuration of the mechanism (zi = 1 or zi = —1). Therefore, the solution 
must come from the last factor and should be valid for all possible values of Z2 . 
Therefore, the coefficients of z 2 ., Z2 and the constant term must vanish, i.e.: 

ai+Mi=0 (33) 

b 1 + b 2 + c + vi + v 2 + w = (34) 

u 2 + a 2 = (35) 

Equation (34) can be written in terms of the design parameters in the following 
form: 

Ji + l 2 2 m 3 + J 2 = (36) 
This solution is physically not possible. 

4.3.4 Reducible case III, Mode A 

This case is completely symmetric with case IV mode A. 

4.3.5 Reducible case III, Mode B 

This case corresponds to the deltoid case with z\ = 1. By replacing z\ = 1 in 
the equation of the angular momentum, K can be written as: 

K = (z 2 - l)(z 2 + l)((ui + v 2 )z 2 2 + wz 2 +v 1 +u 2 ) = (37) 



l 2 Z 



Using the same arguments as in case HB, the term (u\ + v 2 )z 2 + wz 2 + v\ + 
u 2 = should vanish for every unit complex number z 2 , therefore: 

mi + v 2 = (38) 
w = (39) 
v 1 +u 2 = (40) 
Equation (39) corresponds to l\J^l\J 2 — which is physically not possible. 

4.4 Reducible case IV, mode A 

Here we have l 2 — d and l\ = 1% =: /. This is the second case where we get 
solutions that are physically realizable. Again, we introduce the parameters q% 
and Ji and eliminate Ii,ri,pi for i = 1,2, 3. Here are the balancing conditions. 

qi = q z - lm 3 , q2 = ~J<l3, (41) 
I 2 - d 2 

J 1 = J 3 -l 2 m 3 ,J 2 = q 3 -J 3 . (42) 

It follows that qi, q 2 , q 3 must be real. The parameters fulfill again the inequality 
constraints 

m t > 0, J l m l - | ft | 2 > (43) 
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Kinematic constraints 


Kinematic mode 


Static balancing 


Dynamic balancing 


l\ = I2, h — d 


Zl ^ z 2 


Fi = F 2 = 


possible iff d > y/2 Z 2 




Zi = z 2 


Fi +F 2 = 


no 


h = d, I2 = I3 


Zl ^1 


Fi = F 2 = 


possible iff d > \[2 I3 




zi = 1 


F 2 =0 


no 


h = d,l\ = I3 


Z2 + -1 


Fi = F 2 = 


possible iff d > V2 l 3 




z 2 = -1 


Fi = 


no 


Otherwise 




Fi = F 2 = 


no 



Table 6: Complete balancing constraints for planar four-bar mechanisms. 



for i = 1,2,3; again it follows that J\ and J 2 must be positive. By equation 
(42), we get an upper bound for 7713, which must be larger than the lower bound 
from (43) (i = 3). This yields 

(J 3 -lq 3 )(J 3 + lq 3 ) >0. (44) 
This is equivalent to the statement J3 > | Z ^3 1 . From (42), we obtain an upper 



bound for J3 , namely 



<r-r 

<?3 



The lower bound must be larger than the upper 



l 

bound, hence (73 < and d > \/2l. 

Conversely, if d > y/2l, then we can choose 53 < arbitrarily and J3 between 
-(lq 3 ) and g — [ i -q3- Then we choose m.3 between the upper and lower bound for 
77J3 derived above. Next, (42) determines J\ and J 2 , which will then be positive. 
Then (41) determines q\ and g 2 , and finally mi and m 2 can be chosen so that 
inequality (43) is fulfilled. 



4.4.1 Reducible case IV, Mode B 

This case is completely symmetric with case III, mode B. Therefore, there are 
no solutions in this case. 



4.5 Reducible case V 

These cases are similar to II-B, III-B and IV-B. Using the same approach as 
above, the mechanism cannot be dynamically balanced. 



5 Conclusion 

The complete charaterization of dynamically balanced planar four-bar mech- 
anisms is given in table 6. Note that these simple balanced mechanisms can 
be combined to build more complex planar and spatial balanced mechanism as 
shown in [9] . 
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